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ABSTRACT 


An analysis of the unsteady aerodjiiamics of bodies with concave nose geometries 
has been performed. The results show that the experimentally observed pulsating 
flow on spiked bodies and in forward facing cavities can be described by the developed 
simple mathematic model of the phenomenon. 

Static experimental data is used as a basis for determination of the oscillatory 
frequency of spike-induced flow pulsations. The agreement between predicted and 
measured reduced frequencies is in general very good. The spiked-body mathematical 
model is extended to describe also the pulsations observed in forward facing cavities 
and it is shown that not only can the frequency be predicted but also the pressure time 
history of the pulsation phenomenon can be described with the accuracy needed to 
predict the experimentally observed time average effects. This implies that it should 
be possible to determine analjiically what the impact of the flow pulsation is on the 
structural integrity of the nozzles on the jettisoned empty SRM-shells. 
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Section 1 
INTRODUCTION 


According to present plans the Solid Rocket Motors (SRM's) of the space shuttle 
launch vehicle are to be jettisoned after burnout and recovered for renewed usage 
(Ref. 1). Consequently, wind tunnel tests have been made of the SRM's in the com- 
plete angle of attack range 0 < a < 180° (Ref. 2). Recent tests (Ref. 3) showed 
that the SRM could have a stable trim point when it was "flying L.okward" at super- 
sonic speed, and that discontinuous changes of the aerodynamic characteristics occur 
in that or-range, 140° < a < 180° (Fig. 1). The discontinuous flow change is associ- 
ated with a-hysteresis and pulsating unsteady flow (Fig. 2), all characteristics that 
have been observed on spiked bodies. Consequently, the spiked body analysis per- 
formed earlier, when determining the effect of the escape rocket on the Apollo- 
Satum launch vehicles (Refs. 4-7), has been reexamined. New analytic methods have 
been developed which can predict the frequency at which large scale spike-induced 
oscillatory or pulsatory flow occurs. Ihis analytic tool has been applied, after some 
modification, to obtain an assessment of the frequency at which flow pulsations occur 
on forward feeing cavities, e.g. , the SRM-rocket nozzle. A simple formulation of 
the pressure time history during one pulsation cycle has also been suggested which 
provides time average results that agree well with the experimentally observed dis- 
continuous changes in the static aerodynamic characteristics of the SRM rockets. 
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Section 2 
DISCUSSION 


Kabelitz (Refs. 8 and 9) has performed a thorou^ theoretical and experimental 
investigation of the unsteady flow on biconic concave nose geometries. He shows how, 
when the spike is lengthened, the pulsating flow phenomenon (Fig. 3a) is transformed 
to another unsteady flow geometry (Fig. 3b), which he calls "oscillations" and Album 
(Ref. 10) has referred to as "flexions." l^belitz also discusses a third form of flow 
unsteadiness, "vibrations," which occur when the spike is lengthened further (Fig. 4). 

2.1 Causes of Flow Unsteadiness 

Maull (Ref. 11) explained the pulsating flow occi.rrence in the following simple 
manner. When the external flow U^ can be turned by an attached conical shock, a 
steady flow condition exists (see top half of Fig. 5). Maull computed the largest 
spike length that would allow this attached conic shock for various shoulder radii on 
the flat cylinder face. He found that this critical geometry agreed rather well with 
the experimentally established boundary between stable and pulsatory flow (Fig. 6). 
When the required flow turning angle is larger than what can be accomplished by such 
an attached conical shock, a detached strong shock is formed which causes most of 
the shear layer to be turned back into the recirculatory region (see bottom half of 
Fig. 5). As a result the recirculatory region starts to grow and the pulsating flow 
is established. 

It may seem overly simplistic to use this purely inviscid flow concept to deter- 
mine when a viscous flow phenomenon starts. However, Elfstrom (Ref. 12) has used 
the same concept to determine when incipient separation starts in high Reynolds 
number flow; it starts v/hen the required flow turning angle exceeds that allowing an 
attached oblique or conic shock. Bogdonoff et al (Ref. 13) show that Elfstrom's 
sinqjle inviscid concept agrees better with t oerimental data than other more involved 


2-1 


LOCKHEED MISSILES & SPACE COMPANY. INC. 



I 


i 


LMSC-D057194 


mathematical models (e.g, , Ref, 14). The same simple inviscid flow concept also 
predicts when shock-induced separation will occur on a delta wing at supersonic 
speeds; it occurs when the leading edge shock becomes detached, thereby creating an 
embedded supersonic flow region, that is terminated by a normal shock, which usually 
will be strong enough to cause the boundary layer to separate (Refs. 15-17). 

Kabelitz (Ref. 8) speculates that the reasoning behind Maull's criterion for flow 
unsteadiness was that the flow aft of the reattachment shock(s) (aft of WE in Fig. 6) 
goes subsonic, which of course almost coincides with the shock detachment event 
(Ref. 18), Kabelitz (Ref. 8) has himself performed a more elaborate analysis. The 
stationary closed separated flow domain is assumed known, defined by available ex- 
perimental data. The stability boundaries for this closed separated flow region are 
determined using a linear perturbation analysis in which the "unknowns" are lumped 
into one global damping parameter (X^). It is shown that the boundary between oscil- 
latory and pulsating flow is rather insensitive to the value of this parameter A 
and that the predicted boundary agrees well with that determined experimentally . 
However, the boundary between stable and oscillatory flow cannot be determined be- 
cause it is very much dependent on the unknown value of the damping parameter A^. 

2.2 Determination of Oscillation Period 

Frames b through f in Figure 3a show how the recirculatory region grows as 
a result of the free shear layer being turned back into the recirculatory region by 
the strong (reattachment) shock. When the front of this growing separated flow region 
reaches the spike tip the separation region collapses and is "spilled" downstream over 
the cylinder shoulder in the form of a ring vortex (Ref. 19). The recirculatory flow 
build-up then starts all over again as is shown in frame a of Figire 3a. Thus, the 
time period for one cycle of pulsation is composed as follows. 


The analysis is limited to pointed, concave, bi conic forebodies, 
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is the mean upstream convection velocity of the mass flow added to the 
recirculatory region by the strong reattachment shock interaction. After reaching the 
spike tip the separated flow region is "spilled" downstream with the free stream 
velocity. 


The Strouhal number defined by Eq. (1) is 



1 + 



( 2 ) 


The approach used earlier in the analysis of spike-induced unsteady aerodynamics 
(Refs. 4-6) will be used to determine the convection velocity U^. Thus, how 
varies with spike length is going to be determined through the spike efficiency as a 
forebody drag reducer. 


The stable and imstable flow conditions for a critical spike geometry (the crosses 

and circles in Fig, 6) are illustrated in the top and bottom halves, respectively, of 

Figure 5. Thus, the overflow mass flow being turned back in the unstable case is 

transported upstream with the steady state convection velocity, existing at an 

earlier fme instant. It is assumed that the velocity in Eq, (2) is proportional 

to the mean convection velocity, U , for the corresponding fictitious steady flow 

uc 

condition . This steady state convection velocity U will decrease with decreasing 

uc 

spike length since the back-flow-driving pressure rise uirough the reattachment shock 
is decreasing (Fig, 8), In Figure 9 the hypothetical steady state flow case is repre- 
sented in idealized form. It is assumed that for the spike lengths of interest the flow 
profiles in the top half of Figure 9 remain similar. An equivalent flow geometry, 
from drag reduction standpoint, is sketched in the bottom half of Figure 9. Using it 
the maximum drag reduction, obtained with the longest spike that will retain flow 
separation from the spike tip, can be expressed as follows; 

(ac. \ ~ (d U ^ (3) 

\ \/MAX \ °/MAX 


If: 

iilbelitz (Ref. 8) lias shown that application of steady flow parameters to the unsteady 
flow works well. 
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Neglecting density variations between different flow profiles the drag reduction 
for a less than optimum spike length can be written 


AC 


~ (^o) 


^ \2 ^ 
MAX 


(4) 


This is due to the fact that the flow reattachment on the cylinder face remains at 
the same near-shoulder locatiou also for the shorter spike. 


According to the assumption of similar flow profiles the following holds for real 
and idealized flow geometries . 


(«uc)MAx/(“e)MAX = '^uc^e 

{\c)mAkK ' 

Eqs. (3) and (5) give the following definition of U/U^ 

il/2 


U/U = 

oo 




'MAX 


With Uy ~ Eqs. (5) and (6) give 


U 


u/(“u)mAX ' 


"\/K) 


MAX 


1/2 


Finally, combining Eqs. (2) and (7) gives 


«LS 

f ^<=Ao 


1/2 1 

!■ 

(®Ls)oPT 

[ ('^^ao)maxJ 



. (^u)mAX 



u 

u 


^ u 

oo 

X 

(^u)mAX 


(5) 


( 6 ) 


(7) 


( 8 ) 
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Di order to obtain information about ^^^q/(^^Ao)mAX ^ study 

geometries that allow steady flow for the (Lg/D) - ^ange of interest^. Spikad bodies 
with rounded (e.g. , spherical) shoulders form one such family (Ref. 10). From an 
earlier study (Refs. 5-7) the results shown in Figures 10 and 11 were obtained. 
Using the linear representation in Figure 11 one can write Eq. (8) as follows: 


^LS* ^ (^Ls)oPT ^ 

®LS ®LS 

For flat-faced bodies v/ith 
Pulsating flow data for such bodies (Refs. 20 and 21) give the results shown in Fig- 
ure 12. The data have been normalized to the reattachment conditions for the optimum 
spike length using "q. (9). The data are seen to scatter around a straight line. Tne 
decrease of Sj^g* with increasing (Lg/D ■ ^^jj/D) reflects the fact that the second 
term within the "arrow-bracket" in Eq. (8) decreases as U^/U^ MAX ^'^*^^^ 2 .ses with 
increasing spike length. The results in Figure 12 together with Eq. (9) define Sj^g 
as follows: 


kJl-T) 


il/2 


(9) 


thin pointed spikes L/D can be substituted with L^/D. 


S 


LS 


0.32 




TTT 

OPT 




1/2 r 


1-0.14 



( 10 ) 


The experimental results in Figure 10 can be approximated as follows: 



1.6 (1 + 0.1 M ) : < 11.3 

' oo' oO 

3.4 : M >11.3 


( 11 ) 


For the bluff geometries producing pulsating flow a good approximation of the 
(no spike) shock detachment distance is the following (Ref. 22). 



where p/p is the density ratic through a normal shock. 

OO 

If: 

Kabolitz (Ref, 8) lias shown f'at application of steady flow parameters to the 
unsteady flow works well. 


( 12 ) 
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Combining Eqs. (10) - (12) gives the Sj^g - predictions shown in Figure 13 for 
M =2, 4, 6, 3, 10, and oo. The agreement with avi i 'able experimental data 

OO 

(Refs. 11, 20, 21, 23-26) is generally very good. It is oal> Ivlaull's single data point 
(Ref. 11) that shows any substantial deviation, llie reason for this is not known. The 
results in Figure 13 are all for pointed thin spikes. If the spike has a blunt tip, e, g. , 
a planar cut, the measured frequencies (Ref. 24) are higher than what is predicted by 
Eqs. (10) - (12) . The reason for this is that the reactachmant conditions are those 
corresponding to a longer pointed spike (see sketches ir Fig. 14). In this case, and 
also when a "hammerheaded” spike tip is used, one needs forebody drag measurements 
in lieu of Eqs. (10 and (11) for the determination of S. _. 


It is deal that the similar profile assumption used to deri/e Eq. (10) is not valid 
for short spike lengths . When Lg/D exceeds the established reeirculatory 

flow region starts out with a finite upstream convection velocit\- and not with = 0, 
which is the value used in Eq. (10). Tha.i io, the predicted Sj^g - values in Figure 13 
should be finite, not zei ), when (Lg/D - approaches zero. One way to accom- 

plish this is to extend the predictions tangentially below a certain shock penetration. 


Judging by the ex^vcrimental data in Figure 13 it appears reasopable to make this 
linear extension fer (Lg/I^ ~^sh^^^ “ 0.25. Thus Eq. (10) is modified as follows: 


0.32 


’LS 


3PT 



'/2 


/l-c a . 

IS s.v 

1 n ^sh\ 

1 - 0.14 y - p j 


\D ■ rr 



sh 


0.26 


1 + 3.46 


V D D / 


J 


D 

D 


D 

A 


> 0.25 


(13) 


sh 


D 


0.25 


The similar flow profile assumption also breaxs dov/n at the other end of the 
(l/,/D - ^gj^/D) - range. For deep spike penetration the reattachment conditions 
approach a constant value. F,q. (2) shows tliat for = con..tant S, g = constant. 


jj. 

It should be noted that J^Iaull (Ref. 11) used a pointt^d spike, and that no corner 
roundness is present for the daUi point shown in Fig. 13, 
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This constant value for large spike lengths is Sj^ . Figure 11 shows that in the 
liner r approximation = (^^Ao)mAX 



(14) 


Together with Eqs. (10) and (11) Eq. (14) gives 

Sj^* > 0.26 - 0.024(^j|^^^ (15) 

Combining Eqs. (11) and (15) gives Sj^g as a ivJiction of Mach number (Fig. 15). 
As this Sj^ is for flat-faced cylinders. Figure 15 represents the maximum 
Strouhal number for all spike-induced flow pulsations. Figure 11 shows that the de- 
viation from the linear relationship, and the associated similar flow assumption, 

Eq. (10), begins at 



(16) 


Letting Sj^g reach Sj^g in a smooth high order fusion for Lg/D - = 

(Lg/D)^p^ gives the following modification of Eq. (13) 



i 
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0 . 059(0 )o PT \) 

0.45 (Lg/Djop^ 

0.260 - 0.0535 
0.214 - 0.030 

Figure 16 shows the predictions by Eq. (17) compared with the experimental data 
'Refs. 11, 20, 21, 23-26). As most of the experimental data fall in the region for 
which Eq. (10) is valid, 0.25 < (Lg/D - improvement 

in agreement t*^,tween experiment and theory is not very striking when comparing 
Figs. 13 and 16. If the main body face is not flat but of some convex geometry, an 
effective spike length L can be defined as the axial distance between spike tip and 
reattachment (see inset in Fig. 10). hi this case the reattachment diameter D may 
also be significantly different from D. Thus L/D has to be used instead of Lg/D. 

2.3 Effect of Forward Facing Cavity 

When the face of the inain (spike-supporting) body consists of a thin-lipped forward 
facing cavity, the pulsating flow condition persists even for long spikes that would 
produce stable flow on a regular flat faced body (see Ref. 27 and Fig. 17). It is clear 
that the depth (h^) of the cylindrical cavity should be included in the characteristic 
length when defining a Strouhal number for cavity-spike geometries. This effective 
Stioaiial number (Sj ^ should approach that for regular spiked bodies (Sj^) when the 
caviU depth (h^) goes to zero, n.b. provided that pulsating flow would exist for the 
regular spiked-body geometry, hi view of this it appears reasonable to let this union 
between Sj^ and take place at the break-away point defined by Eq. (16). 

The break-away for the regular spiked body occurred because Sj^ was approach- 
ing an upper limit (Sj^g ) when the spike length was increased further, a limit that 
never is realized because the flow becomes stable before the optimum spike length is 
approached. For the cavity-spike geometry the break-away will be in the opposite 
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n = 


X. = 


’LS 
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direction. As more and more mass flow is captured by the sharp cavity rim with 

— # 

increasing spike length, in £q. (2) continues to increase . If one assumes that 
this increase occurs in a (locally) linear fasion the break-away takes the simple 
tangential form and Eq. (13) is modified as follows: 



Figure 17 shows that the predictions through Eq. (18) are in good agreement with 
the experimental results obtained by Hermach et al (Ref. 27). Also shown in Fig- 
ure 17 is one data point for the empty SRM shell "flying backwards." hi what follows 
it will be explained how this data point was introduced into the graph of Figure 17. 

Johnson (Ref. 28) has performed tests at high Mach numbers, > 21.6, of an 
unspiked forward facing cylindrical cavity. He shows that the flow is stable and the 
shock configuration can not be distinguished from that of a flat faced cylinder. This 
is for O' = 0 and no mass addition into the cavity . When mass addition is provided 
through an orifice in the center of the cavity bottom, this "fluid spike" starts a 
pulsating flow phenomenon. In difference to that observed by Hermach et al on the 

tf. 

‘ Provided that the cavity depth is fully effective, i.e., for h^/D < 1. 
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spiked cavity (Ref. 27), this pulsating flow phenomenon is asymmetric and nonpeririic. 
Johnson describes it as follows: "shock bulges out on one side but the asymmetry 
rotates. " The statement about non-periodicity should be judged against the fcict that 
only high-speed motion pictures were used to arrive to this conclusion. Demetriades 
(Ref. 26) has observed a similar asymmetric pulsating flow phenomenon on a stubby, 
tension-shell -like, spiked body. He finds, using fluctuating pressure data, that also 
this rota ting-pulsating flow phenomenon definitely is periodic. 

When Johnson (Ref. 28) changed the cavity geometry to a conical shape, the flow 
again was steady without mass addition but became unsteady when the "fluid spike" 
was present. In this case, however, the asymmetric flow phenomenon does not 
rotate. Johnson's description is: "part of the shock is fixed, extending into the 
cavity; other part of shock snaps in and out." These same words also describe what 
one sees in the flow photographs for the pulsating flow phenomenon on the SRM nozzle 
in Figure 2, and on the hollow skirt of a Polaris re-entry body in Figure 18 (Ref. 29). 
In these latter cases the asj-mmetric flow pulsation is started when at a certain flow 
inclination the windward branch of the bow shock jumps inside the conic ca\ 'ty (see 
Figs. 2 and 18). The sketch in Figure 19 illustrates how this flow pulsation at non- 
zero angle of attack (Fig. 19a) is the "asjTnmetric equivalent" to the pulsating flow 
existing at a = 0 for a spiked cavity (Fig, 19b). The sketches show the hypothetical 
unstable reattachment conditions starting the pulsating flow. 

hi the case of the spiked cavity the most forward extent of the pulsating separated 
flow region is determined by the spike tip. This is als.. true for a regular spiked 
body, in which case the smallest forward extent is obtained when the spike just pene- 
trates the no-spike bow shock, i.e. , Lg > for pulsating flow. Figure 2 shows 
that the distance between the most forward and aft positions of the pulsating strong 
shock (with associated separated flow region) is approximately equal to the shock 
detachment distance for stable flow at a = 0. It therefore seems reasonable to as- 
sume that the "effec ive spike length" is Lg = k A^j^. The results in Figure 2 give 
k = 1, and for lack of more extensive test results this value will be used in the present 
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analysis . * Thus , the SRM-nozzle data point in Figure 17 was obtained fi'om the 
measuJ’ed Strouhal number = 0.18, based on the SRM^diameter, by convert- 

ing i1, hi an effective spike Strouhal number throu^ the following relationship. 

®LS “ '^d'sRM 

## 

^sh^ ‘ ^ is obtained from Eq. (12). Dj is the mean internal nozzle diameter and 
Dg is the maximum external body diameter in the exit plane (i.e. , including skirt). 

.4 Control of Pulsating Flow 

Ic may be desirable to change the frequency of the flow pulsations to avoid reso- 
nance with a structural frequency. How this can be done by changing spike and/or 
cavity geometry is shown by the analytic formulas developed earlier in subsections 2.2 
ard 2.3. If the change is large enough the more dangerous pulsating flow phenomenon 
tx insfcrms into oscillations or possible to steady flow* (see Fig. 4). Rounding the 
face shoulder (Ref. 11) or decreasing the cone angle from = 90° (Ref. 8) are 
ei.ective means for stopping the flow pulsations. When these straightforward control 
m as'ures do not suffice, one has to use more sophisticated means. Obviously the 
reattach ment pro- ^.ss and associated unstable flow phenomena are very much dependent 
upon I'le character of the impinging free shear layer. If the shear layer profile is 
stretched, the re attachment conditions at the shoulder of a spiked body may approach 
those Ihr a small protuberance in a thick boundary layer (Refs. 30 - 32). Such a 
change of the free shear layer profile is obtained by increasing the Mach number and/ 
or decrer.sing the air dersity. 


*Alt*' jugh k should be dependent on at what a' the pulsations occur, the a '-range 
F^ems to be limited to 35° < a'< 45° (Figs. 1 and 2), and a k-value independent of 
’ car be used. 

Volume wise, t!^t is. It is the diameter marked D in Fig. 17. 

*Kabelitz' vil rations are included in the "steady flow" nomenclature used in the 
prese ' i port. 
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Figure 20 shows how the reattachment configuration on a flow separation disk 
undergoes this type of shear layer thickening when the Mach number is increased 
above M„ = 5 (Ref. 5). The absence of the strong reattachment shocks at the disk 


00 


rim for M >5 rignifies a loss of the discontinuous or hi^y nonlinear disk loading 

which causes dynamic instability at < 5. Reeves (Ref. has determined the 

critical shear layer thickness for this change from stable to pulaating spike-induced 

separated flow using the experimental results obtained by Maull, Mair, and Wood 

(Refs. 11, 23, and 34 respectively). (See Fig. 21). The results obtained by Maull 

and Wood are for laminar flow, and for Ljg/D « 2.3. Thus, in Figure 21 is 

in that case a measure of the voiiScity of the lamixiar shear layer at reattachment. 

If one considers the spiked-body flow as a special case of cavity flow, one would 

expect the shear layer to have a tjvo-diraensional character. That is, it should in- 
1/2 

crease as x , just as the laminar boundary layer. The laminar boundary layer 
thickness on a cone varies with edge Mach number as follows (Ref. 35). 


6{M^) = ^(0) (l+ K^M^^ 


( 20 ) 


which suggests that 

^SL <^e> 
D 


^SL 

D 


1 + 




M^ 


( 21 ) 


For the solid line curve shown in Figure 21 3 gj^ (0)/D = 0.0165 and Kgj^ = 0.11. 
Thus, Figure 21 indicates that for L„/D = 2.3 and 0^ = 90° the critical value of the 
she? ' Isy^r vorticity is ind»_ pendent of free stream Mach number. The corresponding 
critical ahear layer edge Reynolds number at reattachment is 


(•^sUcrit = 0-29 *»•»«) '■If' 


It is shown by Kabelitz (Ref. 8) that for Mach numbers at or below M^ = 2 the 
experimental data scatter as much as a factor of 2 or more in regard to the value of 


The static effect is the opposite, i.e., statically stabilizing, which is explained at 
length in Reference 5. 

Lg/D is an important parameter in addition to its influence on M^ in Eq. (21). 
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(Lg/D)^^^ for similar Reynolds numbers. Thus, it probably is coincidental that 
Mair's data point ,(Lg/D)^^.^ = 1.3 for = 1.96, fits on the curve in Figure 21. 
Another way to accomplish this stretching of the shear layer profile is to use a 
"hammerhead" spike tip and let the mixing with the "dead air" base flow region do 
the work. Beastall and Tlimer (Ref. 36) have demonstrated how this can eliminate 
the flow unsteadiness (Fig. 22). Reeves accomplished the same thing by letting a 
ring suspended in the shear layer (upstream of reattachment) do the mixing (Ref. 33). 

hi the case cf the forward facing cavity, e.g. , the empty SRM-nozzle (Fig. 2), 
control of flow pulsation can also be carried out in a manner similar to what has been 
discussed for the spiked body. Although the frequency dependence on geometry is 
less well established it should often be possible to make sure that resonance with a 
structural mode is avoided by using Eqs. (18) and (19). hi regard to the possibility 
to inhibit or moderate flow pulsations both external and internal nozzle geometry are 
important. It is clear that steady flow can only exist as long as the bow shock stays 
external to the nozzle. The maximum flow inclination for which tliis will occur can 
be increased by "bluntening up" the nozzle exterior. When the unsteadiness finally 
occurs at a higher flow inclination than before, the associated unsteady loads should 
be less than before. Using the projected free stream velocity component on the noz- 
zle axis as a measure one would expect the sound pressure level to vary as the n'th 
power of the cosine of the flow inclination, where n > 1. 

The internal nozzle geometry has also a large influence on the pulsating flow 
phenomenon. It was discussed earlier what the effect was of going from a cylindrical 
to a conic cavity. Johnson (Ref. 28) also showed that a hemispherical internal nozzle 
shape caused flow unsteadiness even without mass addition. The unsteadiness was of 
the rotating type observed on the cylindrical cavity with mass addition except for 
"snaps into stable position for short intervals." With mass addition in the center of 
the cavity bottom a seemingly steady asymmetric bow shock configuration was obtained 
(at a = 0) which looked very similar to the "end-stroke position (most forward) of the 
bow shock geometry obtained on the SRM-nozzle at large flow inclination (Fig. 2). 
Johnson's results are for extremely high Mach numbers, M^ > 21. However, simi- 
lar results for hemisphere cavities have been obtained in other investigations at 
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supersonic Mach numbers (Refs. 37 and 38). At = 1.98 Ck)oper et al (Ref. 37) 
found a pseudo-stable flow condition, i.e. , the bow shock would stay in the symmetric 
or asymmetric configuration for long time periods with the change between the two 
states coming at random intervals (Fig. 23). Markley (Ref. 39) could, however, 
only obtain the symmetric stable shock configuration. At = 4.95 Cooper et al 
obtained the oscillatory asymmetric pulsating shock configuration that Johnson ob- 
served (Ref. 28). The Strouhal number based on face diameter was Sjj = 0.099. 


Demetriades (Ref. 26) and Cassanto et al (Ref. 40) have tested tension shell type 
geometries tliat combine some of the geometric features of spiked bodies and forward 
facing cavities (Fig. 24). The data points for the regular symmetric pulsating flow 
phenomenon agree weU with the predictions from Eq. (17), as does also the experi- 
mental data obtained by Baltakis (Ref. 41). Both Demetriades and Cassanto et al 
observed also an asymmetric pulsating flow phenomenon very similar to that discussed 
here earlier for forward facing cavities. The frequency of this rotating (according to 
Demetriades) or rocking (according to Cassanto et al) pulsating flow phenomenon was 
observed to be half that of the S 3 onmetric pulsating flow. Using Eq. (18) one computes 
for a hemispherical cavity at =4.95 the Strouhal number Sj^g =0.088. Using 
half this frequency^ gives Sjj = 0.088 which is in rather good agreement with the 
value Sp = 0.099 measured by Cooper et al (Ref. 37). 


Stallings and Burbank (Ref. 38) observed at =2.49 and =3.5 the same 
pseudo-steady bow shock condition that Cooper et al (Ref. 37) found at =1.98, 
and obtained the asymmetric flow pulsation at =4.44 for the whole angie-of- 
attack range tested, 0 < a <15°. However, they did not measure the frequency of 
the flow pulsations. Tliis unsteady flow phenomenon is of great concern when design- 
ing a drag device for use at supersonic speeds (Ref. 42). In that case the design 
solution is to combine the "hammerheaded" spike tip effect (Fig. 22) with an opening 
for partial through- flow at the bottom of the cavity. Although it is obvious that viscous 
effects play a dominant role in setting up these unsteady flows, there is some indica- 
tion that unstable inviscid flow processes also play a role (Ref. 43). 




Noting that 





I in Eq. (18). 
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In the case of the aiM-nozzle "flying backwards" (Figs. 1 and 2) some "throng 
flow" at the base of the cavity is a possibility. However, it will more than likely be 
of much more complicated character than in the case of the "through flow" parachute 
(Ref. 42), From geometric conside .'ations one is likely to incorporate a new unsteady 
phenomenon, that associated with the so-called resonance tube (Ref, 44). This will 
also be the case if the cavity is not shallow, as has been assumed in the discussion so 
fe,r. Figure 25a shows a series of flow photographs illustrating the various stages in 
one cycle of the oscillation taking place in the two-dimensitmal resonance tube used by 
Vrebavolich (Ref. 44). The "flow trip" starting the oscillation, or at least greatly 
amplifying it, is a two-dimensional spike or what Vrebavolich calls a "wing-trip." 
Frame 1 shows the "swallowing" of the normal shock with completely separated flow 
up to the spike tip. In Frame 2 the spike flow is almost fuUy attached while the shock, 
after reflection off the bottom of the tube, is moving upstream. This is the compres- 
sion phase of the oscillation (see Fig. 25b) . hi Frame 3 the shock is being "spilled" 
starting the spike flow separation which in Frame 4 has reached the spike tip and 
remAU'.s there in Frame 5, while the tube is still "spilling" air. hr* ^rame 6 the ex- 
pansion or "swallowing" phase is just about to begin. Figure 25b shows that the strong 
shock reached the spike tip between Frames 5 and 6, as is indicated by the steepening 
of the spike tip shock (strength) (see top graph in Fig. 25b). 


Comparing Figure 25a with Figure 3a one can see great similarities. The big 
difference is that in Figure 25a, i.e. , in the presence of a very deep forward facing 
cavity, tiie period of the oscillation cycle is determined by the filling and spilling 
process of the tube or cavity. This time period is that corresponding to the organ 
pipe oscillation with the wave length equal to four tube lengths . Thus, the core- 
sponding Strouhal number is in this case 



U 

00 



1 + 


(y - 1) 


1/2 


( 22 ) 


The numbers in Fig. 25b indicate where the flow pictures in Fig. 25a belong. 

## 

This is with the tube closed at the downstream end. With it open the wave length for 
the ground tone would be two tube lengths. 
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With h- = 10 inches the measured frequency (Ref. 44) gave f h^/U = 0,122 
whereas Eq. (22) g..es f h^/U^ = 0.130. Using Eq. (22) one obtains for Sj^g 



hi Figure 26 the experimental results for M =3.5 and L_/D = 0.9 (Ref. 26) 
are compared with the predictions obtained from Eq. (18) and Eq. (23). It is obvious 
that for shallow cavities, h^/D < 1, the convective processes dominating the spiked 
body flow pulsation still are in control. For deep cavities, h^/D > 1, the acoustic 
fill-spill cycle of the cavity, the organ-pipe phenomenon, takes over the control. It 
is worthwhile noting that the asymptotic value Sj^ = 0.121 is the constant value that 
would apply if one assumed that the flow recirculation on the spike also was acoustically 
controlled. It is clear that this cannot be the case, especially when considering the 
rest of the experimental data (Fig. 17). One would be inclined to conclude from Fig- 
ure 26 that as long as the convective flow processes on the spike anu cavity inlet 
(shallow cavity) are slower than the acoustic processes in the rest of the cavity (deep 
cavity), the flow oscillation frequency is given by Eq. (18). Conversely, when the 
acoustic cavity process takes the longest time (h^/D >1.15 in Figure 26) the oscil- 
lation frequency is given by Eq. (23) or Eq. (22). 

Figure 25a illustrates the situation often existing for supersonic inlets , a simi- 
larity that has been utilized to analyze inlet buzz (Ref. 45). To wnat extent it applies 
to the interference flow field between or biter and booster (HO tanks and SRM's) on the 
space shuttle launch vehicle is a moot question because of the spanwise spillage and 
the presence of exhaust plumes (Fig, 26). The flow in deep forward facing cavities 
can also become oscillatory without using the spike as a tripping device for the un- 
steadiness. The higher than stagnation value heating of the bottom of these resonance 
tubes have very useful practical applications, and extensive studies have been made 
of them. A recent review of the field is contained in Reference 46, which also shows 
that the oscillatory flow exists even at high subsonic speeds, i.e. , in absence of the 
strong moving shocks hi !he bibe which have been thought to be necessary. The 
authors have no explan, ition for this but speculate that there must be some interaction 
between viscous and acoustic processes (also in absence of shocks) that sets up the 
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oscillation and associated heating amplification, which in their test was as high for 
this subsonic condition as for the various supersonic cases investigated, hi the un- 
steady flow occurring in sideward facing cavities the situation is somewhat similar to 
that for spiked forward facing cavities, i.e. , convective shear layer processes compete 
with cavity acoustics in regard to the control of (he oscillation and its frequency. 

(See Ref. 47 for a thorough and up-to-date review of the field.) TTie flow oscillations 
are only observed in what Charwat et al (Ref. 48) have termed "open cavities," i.e, , 
cavities in which flow reattachment occurs only on the downstream cavity shoulder . 
Thus, it is not surprising to find that a softening of this reattachment through slanting 
the aft cavity wall back; -ards is the most effective means of suppressing the oscilla- 
tions (Ref. 47). The same was true for the spike induced oscillations (Refs. 8 and 11). 


1 


This sets an upper limit on the length-to-depth ratio for the cavity. 


2-17 


LOCKHEED MISSILES & SPACE COMPANY. INC. 



LMSC-D057194 


Section 3 

UNSTEADY AERODYNAMICS 


earlier studies the effect of the spike-induced separated flow on the aerodynamic 
damping has been determined (Refs. 5, 49, and 50), and in a more recent investigation 
also the forcing function off "hammerheaded" spikes has been studied (Refs. 51 and 52). 
In the present analysis only the unsteady aerodjTiamics of the "backward flying" SRM’s 
will be considered. 

3.1 Cavity Flow Pulsations 

Figure 28 shows the innermost shock position before the start of a pumping cycle, 
hi Figure 28a the real flow is sketched with internal oblique shocks non-parallel with 
the nozzle walls and with the normal shock a finite distance above the nozzle bottom, 
hi Figure 28b the idealized flow is shown with oblique shocks parallel with the walls 
and the normal shock an infinitesimal distance above the "floor." This is the 
Newtonian type idealization. Using the definitions in Figure 29 and the generalized 
embedded Newtonian theory of Reference 53 the analysis can proceed in the following 
manner. 

Assuming a two-dimensional Upe behavior for the internal flow one can according 

2 2 

to Reference 54 assume that p U = p U and that C is a function of x only. 

oo <30 po 

Thus, with Vj^/U defined as before (see Ref. 53 with ~ 
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When the strong shock has passed upstream over the area element the pres s 
given by the normal shock relations (Ref. 18). hi the case of steady hypersonic flow 
the shock strength is 

% ‘ rU ^ rh ( 25 ) 


In the nonsteady case the shock strength can be written 



where > 0 for the pumping phase when the shock moves upstream and 

AMjjg < 0 for the expansion phase when the shock moves downstream. For deeper 
cavities the results obtained by Vrebavolich (Fig. 25 and Ref. 44) indicate the follow- 
ing values Ca^ and Ca^ for the shock strength during the upstream and downstream 
movement respectively. 


0 < Ca , < 0.25 C 

pd ^ po 


CAj ^ Ca Ca 
pd pu po 


(27) 
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Tlius, the pressure on the internal nozzle wall will vary with time as is illustrated 
in Figure 30. to the steady case the pressure is Ca^ tlu:oughout the nozzle, where 
< 0^^. Hie pressure change in the unsteady case fi-o.. .ae steady state level is 
AC^ = where is given by either of Eos. (21), (26), or (27) depending 

upon where m ttie oscillation cycle the nozzle is (Fig. 30). The corresponding change 
in the strip loading is (see Fig. 29). 


dAC. 
A 

dr 


dAC^ 

dx 


2r 

— X 


j 




' °go> ■ ‘i - * < ‘2 
(Cv-(>)d# 1 ‘2 =*<‘3 


(Ca - Ca ) sin <t> d(b ; t<t<L 
' pu po' ^ ’ o 1 


(Ca^ - Ca^) sin (t> 6 (p ; ^ t < t^ 


^^p ~ ^po^ sin 0 d <^) ; t^ < t < Ig 


(28) 


(29) 
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Integrating Eqs. (2^) and (29) utilizing Eqs. (24) - (27) gives 


Ca - Ca„ 
po pu 

Ca - Ca , 
po pd 






^ |(f * «”) R =*”' «R 

12 2 \ 1 

+ -; 7 sin a' CCS 0^J +-^sin2o;' sin2 0j^cos 0* 


\‘2^ 


t < t. 


-^C x{ 

s Pmax 


; t. ^t<t2 


+ 0* ~ sin 2 0'j ^ sin 2 a’ sin 2 6j 

[ 2 2 
cos a' sin 0^ 

+ sin^ a' cos^ ^') j ’ * 


i' sin^ 0, 


From Figure 30 one can define the time duration during one cycle for the differ- 
ent pressure levels*. 


(tj - y A = (x - (1 + U/u^) 

<‘2 - = (’' * ^si>'*‘c * ^h> 

<‘3 - ‘2’^ ' "'c ■ '‘•^<^0 * *sh> 


Note that the assumption of non-finite distance between normal shock and nozzle 
bottom at the start of the cycle will have a negligible influence on the computed time 
average loads. 
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H the nozzle walls are thin or the nozzle lips are relatively slnui> one can neglect 
(Be&. 53 and 54). This simplifies the computation of the time average load 
change from the steady state oonditicai. Combining Eqs. (25) - (27) with Eqs. (30) 
and (31) gives the following integrated result. 

ttD. 


AC^ = 


TD/ Dioot9„ (i(, i r Afll 

4S ”ll \VJ) 


[if * *’) 

+ sin^ a* cos^ sin 2 a' sin 2 cos | 

4 ^ 6^tan [<C^ Cja»] /(l ♦ ^) {{ P3) 




C 


S (hQ 




0. 


- 2 Sin 2 0’ 


^ - 2 , 2 ^ 
+ sin a’ CJ3 6 


I ^ sin 2 a' 


' sin 2 + cos 0’ 


r 2 , . 

I cos a' SI 


• 2 
sm 6 


R 


^ cos^ <|)'j 


R 


(34) 


The corresponding change of the pitching moment is 


- *CG> 


X = 


2 " 2 ^R 


3 D 
~2 


I 


(35) 
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Eq. (2) defines U/V = U /U as 

^ ' 00 u oo 


U/U 

00 



( 36 ) 


The discontinuous changes of the aerodynamic characteristics caused by the flow 
pulsations are shown in Figure 31, complementing the - results already shown in 
Figures 1 and 2. The magnitude of the - discontinuity is not measurable. How- 
ever the C . - and C - discontinuties are. Figure 32 shows the measured discon- 
tinuities for three different nozzle geometries (Hefs. 2 and 3). Also shown are the 
time average values predicted from the present analysis. The agreement between 
predicted and measured discontinuities is very satisfactory. Althou^ it may be 
desirable to be able to predict these time average effects on the static aerodynamic 
characteristics, the more significant impact of Figure 32 is that it lends credibility 
to the very simple imsteady flow model suggested for the puls, ting flow process . Thus , 
it should be possible to use the present model in an analysis to determine whether or 
not the pulsating flow phenomenon will endanger the structural integrity of the rocket 
nozzle on the jettisoned, empty SRM-sheUs. 


# 


From Eq, (27) the values 



0 and Ca 
pu 



were chosen for simplicity. 
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Section 4 
CONCLUSIONS 


A study of the imsteady aerodynamics of bodies with concave nose geometries has 
shown the following. 

• The spike-induced flow separation remains steadv as long as Uie recompres- 
sion at reattachment can be accomplished by (a series of) oblique shocks. 

• This is the case if either the spike is long enou^ or the surface inclination 
at reattachment is mild enou^. 

• No unsteady spike-induced flow separation exists if the Re 3 molds number is 
below a critical value or the free shear layer is thicknened artificially to a 
corresponding degree . 

• K the fiice is in lined forward, creating a forward facii^ cavity, the unsteady 
spike-induced flow separation persists even for very long spikes. 

• Unsteady flow exists on forward facing cavities even in absence of a spike. 

At zero angle of attack without mass addition it has been observed only for a 
heir''«'pherical cavity. However, at moderately high angles of attack (at 
which the windward branch of the bow shock moves into the cavity) a pulsating 
flow phenomenon exists which is very similar to that observed in cavities 
with a flow separation spike present. 

• Simple analytic means are developed by which the experimeritally observed 
oscillation frequency of the pulsating flow c: n be predicted both for spiked 
bodies and forward facing cavities. 

• The postulated mathematic model for tlie pressure time history of the pulsat- 
ing flow gives time average results that are in good agreement with the meas- 
ured discontinuities in static tests . This strongly implies tiiat the structural 
integrity of the SRM-nozzle in presence of pulsating flow could be determined 
analytically. 
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Appendix A 
NOMENCLATURE 


axial force: coefficient C . = A/(p U ^ /I) S 
speed of sound 

reference length (SRM-diameter) 

body diameter (Fig. 10) 

reattacbment diameter (Fig. 10) 

idealized reattachment parameters (Fig. 9) 

internal nozzle diameters (Fig. 29) 

diameter of flat face area (Fig. 11) 

spike diameter (Fig. 10) 

oscillation frequency 

cavity depth (l.^. 17) 

boundary layer parameter, Eq. (20) 

free shear layer parameter, Eq. (21) 

spike length definitions (Fig. 10) 

spike length for flat faced bodies, Lg = L = L 

Mach number, M = U/a 

2 

pitching moment: coefficient C = M /(p U /2) Sc 

m p oo 00 
2 

normal force: coefficient C„ = N/(p U /2) S 

N *^00 00 

static pressure: coefficient Cp = (p - p )/(p^ ^oo^/2) 
blast wave pressure 
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S 

®D’ ®Ls 
T 



a 

a' 

y 

A 



9 

9 


c 

R 


P 

<t> 


<P' 


Reynolds number 

Re based on D and free stream conditions 
radius of shoulder roundness 
reference area, S = ttc /4 

Strouhal number, S_ = fD/U ,8, = f L„/U 

U 00 J-iS o 00 

oscillation period 
time 
Ume lag 
axial velocity 

velocity component normal to surface element 

parameter defined in Eq. (17) 

axial body-fixed coordinate (Fig. 29) 

angle of attack 

of-complement, a' = tt- a 

ratio of specific heats, 7 = 1.4 in air 

increment 

shock standoff distance 
bcundary layer tiuckness 
thickness of free shear layer 
spike tip cone migle (Fig. 10) 
slope of nozzle wall (Fig. 29) 
air density 
azimuth (Fig. 29) 

0-value for Newtonian shadow, Eq. (25) 
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I 


1 


LM?C-D05Tliv4 


£ 

e 

eff. 

I. i 


MAX 

NS 

OPT 

S 

SL 


SRM 


u 

uc 

0, 1, 2, 3 


<r> 


Sulecripts 

fextemal 

edge of boundary layer or free shear layer 

effective value 

internal 

linear-tangential extension 

maximum 

nonsteady 

optimum 

spike 

shear layer 

Solid Rocket Motor 

upstream 

upstream convection 

phases in oscillation cyde (Fig. 30) 

free stream condition 


(*) 

(-) 

(a) 

(v) 


Superscripts 

value (hypothetical) for optimum spike length 

barred quantities denote integrated mean or time average values 

"roof" indicates value behind a normal shock 

"inverted roof indicates embedded Newtonian value 
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PULSATION 





Figure 4 Classification of Spiked Body Unsteady Flow Fields 
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Figure 6 Critical Spiked Body Geometry 
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Figure 8 Effect of Spike Length 
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Figure 9 Idealized Separated Flow Geometry 
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Figure 10 Optimum Spike Length as a Function of Mach Number 
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Figure 11 Spike-Induced Forebody Drag Reduction 
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Figure 22 Effect of Spike-Hammerheads on ^ike-lhduced Flow Pulsations 
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Figure 24 Flow Pulsations on Probescidean Nose Geometries 


LOCKHEED MISSILES & SPACE COMPANY. INC. 


s 




iteady How ta a Kes<mattce Tube wltt Flow S»puratl<ai 






LMSC-D057194 




% - ‘ 

LMSC-D057194 



limAUiEo mw 


Figure 28 Idealized Sfoving Shock Process in the SRM Nozzles 
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Figure 30 Idealized Pressure Time History 
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Figure 31 Static Characteristics of the SRM with Nozzle at = 3.48 
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Figure 32 Coo^jarisoti Between Predicted and Measured Aerodynamic 
Discontinuities for Three Different Nozzle Geometries 
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